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Abstract 

We study the spectrum of the QCD Dirac operator for two colors with fermions in 
the fundamental representation and for two or more colors with adjoint fermions. For 
Nf flavors, the chiral flavor symmetry of these theories is spontaneously broken according 
to SU(2Nf) — > Sp(2Nf) and SU(Nf) — ► 0(Nf), respectively, rather than the symmetry 
breaking pattern SU(Nf) x SU(Nf) — > SU(Nf) for QCD with three or more colors 
and fundamental fermions. In this paper we study the Dirac spectrum for the first two 
symmetry breaking patterns. Following previous work for the third case we find the Dirac 
spectrum in the domain A <C Aqcd by means of partially quenched chiral perturbation 
theory. In particular, this result allows us to calculate the slope of the Dirac spectrum at 
A = 0. We also show that for A <C 1/L 2 Aqcd (with L the linear size of the system) the 
Dirac spectrum is given by a chiral Random Matrix Theory with the symmetries of the 
Dirac operator. 
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1. Introduction 



Perhaps the most important reason for the current interest in the infrared sector of 
the Euclidean Dirac spectrum is its relation with the chiral condensate by means of the 
Banks-Casher formula Q. In this way the properties of the smallest eigenvalues of the 
Dirac operator may be related to the mechanisms of chiral symmetry breaking and to 
the nature of the chiral phase transition as observed in lattice QCD simulations @, [!|. 
However, this is not the only reason for our interest in this part of the Dirac spectrum. 
What may be more important ultimately is that it is one of the few examples where exact 
nonperturbative analytical results can be derived from QCD. 

One of the main objectives of this paper is to study the QCD Dirac spectrum in the 
domain A <C Aqcd- In particular we calculate the slope of the Dirac spectrum near A = 0. 
This problem was first addressed in a paper by Smilga and Stern They studied the 
slope of the Dirac spectrum for the case of three or more colors with fermions in the 
fundamental representation. In this article we extend previous results j|, |5], || to the 
case of QCD with two colors and fundamental fermions and QCD with adjoint fermions. 
Before discussing our strategy let us first summarize the approach followed by Smilga 
and Stern. They studied the scalar susceptibility which, in terms of the spectral density 
averaged over the Euclidean QCD action, 

P(A) = -**)>> (1) 

k 

is given by 

^=-^r 'W"'-*V ( 2) 

Jo {m z + x z y 

Their idea was to derive K ab from chiral perturbation theory and invert this relation to 
obtain information on the Dirac spectrum. Because of the implicit mass dependence of 
p(A) via the fermion determinant in the statistical average over gauge field configurations 
such inversion is not possible in general. They found that 

N 2 t - 4 

K ab _ s ab ^- — log(m/A), (3) 
N f 

where A is the momentum cutoff of the one-loop integral. Their crucial observation was 
that the contribution to the integral (0) from the chiral logarithms in p(0) cancels because 

p(0)(m 2 -A 2 ) 



(m 2 + A 2 ) 2 



-dX = 0. (4) 



Barring any more exotic A dependence of p(A) such as for example y/mX logm, the loga- 
rithmic mass dependence of the scalar susceptibility results entirely from the linear term 
in p(A). 
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However, there is another approach to extract the Dirac spectrum directly from the 
low energy effective partition function which does not suffer from the above mentioned 
ambiguities in the mass dependence [0, |], IU- The starting point of this approach is the 
so called partially quenched partition function defined as [§, (|, [UJ 



where Jp is the Euclidean Dirac operator and SVm [-<4] is the Euclidean Yang-Mills action. 
The valence quark mass dependence of the chiral condensate is then given by JTTJ, [12j 

**».) = 7<E = \ S,l^log^(m„, J)- (6) 

The spectral density then follows from the discontinuity across the imaginary axis 

p(A)/y = ^(£(iA + e)-£(iA-e)). (7) 
The valence quark mass dependence of the partition function (Rl) can be obtained 



from its corresponding low-energy effective partition function |9], 10, |TJ]] which is 
based on the specific pattern of chiral symmetry breaking. In this way it is possible to 
extract exact analytical information on the QCD Dirac spectrum in the domain near 
zero virtuality. This approach was first followed in |5|, |5| for three or more colors with 
fundamental fermions where in addition to extending the result of Smilga and Stern we 
were able to calculate the spectral density of the Dirac operator on the scale of the smallest 
eigenvalues. In this paper we apply this approach to QCD with two colors with fermions 
in the fundamental representation and QCD with adjoint fermions. 

A much simpler theory with the global symmetries of the QCD partition function 
is chiral Random Matrix Theory (chRMT) 16fl . In this theory the matrix elements 
of the Dirac operator are replaced by Gaussian distributed random variables. It turns 
out that this does not affect the pattern of chiral symmetry breaking. The separation 
into three different symmetry classes can be related naturally to the Dyson index of the 
Dirac operator denoted by the symbol (3. Its value is f3 — 1 for QCD with two colors 
and fermions in the fundamental representation, it is (3 = 2 for QCD with three or more 
colors and fermions in the fundamental representation, and it is (3 = 4 for QCD with 
adjoint fermions. The matrix elements of the Dirac matrix of the corresponding chRMT 
are real, complex and quaternion real, respectively fIB| , |56| . Within chRMT it is possible 
to obtain exact analytical expressions for the spectral density and all spectral correlation 
functions. As we will see in the next paragraph it is possible to establish a domain where 
the spectral correlations of the QCD Dirac operator are given by chRMT. Alternatively, 
the existence of such domain can be understood by means of universality arguments 



17, 18, lg gO, |21|, |2J, gg, |24 
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We can distinguish three important scales in the QCD Dirac spectrum. They are 
most conveniently discussed in terms of the the valence quark mass. The first scale is 
the smallest nonzero eigenvalue of the Dirac operator which is of the order of the average 
level spacing and is given by 

A min ~ (8) 

where So is the chiral condensate in the chiral limit and V = L 4 is the volume of space- 
time. The second scale corresponds to the valence quark mass for which the Compton 
wave length of the associated Goldstone boson is equal to the size of the box f7|, |25], |5|, 

„ J_ (9) 

As was argued in |7], ^| |^, |6| the valence quark mass dependence of the chiral condensate 
is described by chRMT for 

F 2 

m v < m c = — — -. (10) 

The spectral density of the Dirac operator, given by the discontinuity of £(m„), is therefore 
described by chRMT as well in this domain. These results have been confirmed by numer- 
ous lattice QCD simulations ||, 0, H H 13 0. H © H HI HI M> HI H 



instanton liquid simulations (41], |25|, Q and other calculations [43, 44]. For valence quark 
masses below this scale the kinetic term in the effective Lagrangian can be ignored. In 
|| it was shown that in this regime the valence quark mass dependence that follows 
from partially quenched chiral perturbation theory coincides with the result from chRMT. 
The dimensionless ratio of the scales (|J) and fllPp 

m * FV (id 



represents the approximate number of eigenvalues that are described by chRMT. A scale 
analogous to m c occurs in mesoscopic physics where it is known as the Thouless energy 



45, 46 



Based on the scales 

A m in < rn c < Aqcd, (12) 

we thus can distinguish four different domains. With names borrowed from mesoscopic 
physics, these scales separate the quantum domain, the ergodic domain, the diffusive 
domain and the ballistic domain. The valence quark masses can be interpreted in terms 
of a conjugate time. Then h/m c is the time scale for which a particle diffuses over the 
length of the system with diffusion constant given by F 2 /Sq \WA Hj UMl- For more 
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discussion and additional references on impurity scattering in disordered systems and the 
Thouless energy we refer to |^9, 50, 51 . 



In the domain of validity for chRMT it is natural to study the spectral density in units 
of the average level spacing. The corresponding microscopic spectral density is defined by 



15 



ps{u) = ^w p{ w ] - (13) 

In this paper, we will show by means of a perturbative calculation that in the ergodic 
domain the microscopic spectral density that follows from partially quenched chiral per- 
turbation theory coincides with the result obtained from the chiral Random Matrix Theory 
in the same universality class. For smaller values of m v the super-integrations have to be 
performed exactly. This has been achieved |5|, |(J for the case (3 = 2 which is mathemati- 
cally much simpler than the cases (3 = 1 and (3 = 4 H9, Q and complete agreement 
with chRMT was found. It would be interesting to perform a similar calculation for the 
other two well. 

As second important result of this paper is the derivation of the slope of the Dirac 
spectrum for QCD with two colors and fundamental fermions and for QCD with adjoint 
fermions. We will derive this result in two different ways. First, according to the method of 
Smilga and Stern, and second by means of the partially quenched QCD partition function. 
We find that both methods agree thereby excluding a more exotic mass dependence of 
the spectral density which would be allowed by the Smilga-Stern approach. Our result 
for the slope of the Dirac spectrum for Nf massless fermions can be summarized as 

„<(A = 0) = S | W-2)W + fl (14) 
Hy ' 16tt 2 F 4 0N f K ! 

where S is the chiral condensate, F is the pion decay constant and (3 is the Dyson index 
of the Dirac operator. 

In the next two sections we study the symmetry breaking patterns for (3=1 and 
(3 = 4. The slope of the Dirac spectrum for these two cases is calculated from standard 
ChPT in section 4. Partially quenched chiral perturbation theory based on supersymme- 
try is constructed in section 5. Within these effective theories we calculate the valence 
quark mass dependence of the chiral condensate. The computations for two colors with 
fundamental fermions and adjoint fermions are presented separately in sections 6 and 7, 
respectively. Some interesting limiting cases and compact expressions for all (3 are given 
in section 8. In section 9 we show that our results agree with chRMT in the ergodic 
domain. We also find the Dirac spectral density for larger eigenvalues, in the diffusive 
domain, and obtain an expression for its slope. Concluding remarks are made in section 
10. 
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2. Symmetries of the QCD Partition Function with Two Colors and Funda- 
mental Fermions 



The QCD partition function for two colors and fundamental fermions is special in 
the sense that it is invariant under a larger symmetry group than is the case for more 
colors. This well-known fact for QCD with fundamental fermions relies on the pseudo-real 
nature of SUc(2) |54], [55|, |56| , and can be easily generalized to the supersymmetric case 



57] relevant to the partially quenched QCD partition function. The discussion in this 



section is mainly included for pedagogical reasons and part of it was already given in 
reference |58| . 

The fermionic action is obtained by writing the determinant in the partition function 



as a Gaussian integral over the Grassmann fields 
of the 7 matrices with 



7a. 











and (f). Using the chiral representation 



(15) 



and a 



,, - K l,i<Jk) with Ofc the Pauli cr-matrices, the fermionic action can be written as 

Nf 

S F [ 

./-!. 

(16) 




rrif 



aJdu + iA 




For N c 
with the relation a 



2, we have that A 1 :, 



-t 2 A^t 2 (with r 2 = o"2 in color space). Combining this 



OiGuGi we find 



V L d-Udp + iA„ 



^o- 2 r 2 (J M (9 M + iA^a 2 T 2 (j) f L , 



(17) 



where we have used that is ant i- Hermit ian and that the fermion fields are anti- 
commuting Grassmann variables. The fermionic action can thus be rewritten as 



Sp{mf = 0) 




o- M (<9 M + iA^ 






<7 M (<9 M + iA 




Obviously the symmetry group is Gl(2Nf). If the number of left-handed modes is not 
equal to the number of right-handed modes because of the anomaly, an axial U(l) is 
broken explicitly as for three or more colors. 
The mass term is given by 






\m f a 2 r 2 

\m f a 2 T 2 




-^m f a 2 T 2 




0- 2 T 2 (j){ 



(19) 
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Also in this case we expect a maximum spontaneous breaking of chiral symmetry consis- 





tent with the Vafa-Witten theorem [59|. This means that only the subgroup of U(2Nf) 
that leaves both 0_r0_r and invariant remains unbroken. Therefore, only the sub- 

group that leaves 

/ xf \ / n _ _ \ / J./ \ 

(20) 

invariant remains unbroken {4>r4>r is invariant under the same transformations). This 
is the symplectic group Sp(2Nf). The Goldstone manifold is thus given by the coset 
SU(2N f )/Sp(2Nf). 

An analogous argument can be made for bosonic quarks. In that case it is essential 
to allow only symmetry transformations that leave the integrals in the partition function 
(0) convergent. In that case one finds [E7|] that the chiral symmetry is broken according 
to Gl(2,R) -> 0(2). 

The symmetries of the partially quenched partition function for two colors with fun- 
damental matter fields can be obtained along similar lines. One easily convinces oneself 
that the symmetry of the Lagrangian is enhanced to Gl(2Nf + 2|2) from the usual group 
Gl(Nf + 1|1) x Gl(Nf + 1|1) for more than two colors. For convergence reasons the 
symmetry group has to be restricted to a non-compact group Gl(2,R) for the bosonic 
quarks. This symmetry is broken spontaneously by the formation of a chiral conden- 
sate in the boson-boson sector and the fermion-fermion sector exactly as discussed above. 
The remaining symmetry group combines into the ortho-symplectic graded Lie group 
OSp(2Nf + 2 1 2). The Goldstone manifold is therefore given by |]57| the maximum Rie- 
mannian submanifold of the coset space Gl(2Nj + 2|2) /OSp(2Nf + 2|2). This results 
in a integration domain for the the fermion-fermion block given by the symmetric space 
U (2Nf) I Sp(2Nf) and an integration domain for the boson-boson block given by the sym- 
metric space Gl(2,R)/0(2). 

An explicit parameterization of this Goldstone manifold is given by 

V T GV, (21) 

where 



G = 

and V is parameterized as 

V 











JV>+1 














h 



(22) 



^ Uff di «2 

Pi a b | . (23) 
\ P 2 c d 
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Here, oti and $ are Grassmann variables, Uff is an (2Nf + 2) x (2Nf + 2) unitary matrix 
and a, b, c and d are real numbers. 



3. Symmetries of the QCD Partition Function for Adjoint Fermions 

For QCD with fermions in the adjoint representation the Dirac operator is given by 

JP = + f ahc l»A atl , (24) 

where the f abc denote the structure constants of the gauge group. In this case, the operator 
7274 {Jp + m ) is antisymmetric and therefore the square root of the fermion determinant 



can be represented as a Grassmann integral [60 in the following way 



det 1/2 ( 7274 ( J D + m)) = 

¥r \ ri ^ ( rn a,(d, + f°» c jA aiM ) \ ( ^ 



dipFidipL exp 



(25) 



For simplicity we will consider only even Nf. The flavor symmetry group is SU(Nf) which 
acts on the Majorana fields as 



4>r\ u o \( fa 



L 



(26) 



This symmetry is broken by the chiral condensate which, with maximum breaking of 
chiral symmetry consistent with the Vafa-Witten theorem, is given by (the color index is 
denoted by i and the spinor index by a) 

2 Nf 

So= F E(^(4ai4i + ^aAi)>' ( 27 ) 
iV / f=l 

and is invariant under an O(Nf) subgroup of SU(Nf). The Goldstone manifold in this 
case is thus given by SU(Nf)/0(Nf). 

The square root of an anti-symmetric matrix cannot be represented as an integral 
over real bosonic variables. Therefore, the above construction fails for bosonic quarks. 
However, it is possible to represent the inverse determinant as an integral over complex 
variables. If we use that 

0l<3-+(<9 m + f abc ^A a ^ R = (a 2( f> R )&^ + r fcc 7 ^>20L, (28) 

we observe that the flavor symmetry group of the bosonic quarks is doubled. In this case 
it turns out f57j that the chiral symmetry is broken according to U*(2,R) — > Sp(2). 
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The symmetry group of the partially quenched Lagrangian for matter fields in the 
adjoint representation is Gl(Nf + 2|2). Because of convergence reasons only a U*(2,R) 
subgroup of the boson-boson block of Gl(Nf + 2\2) is a symmetry of the partition function. 
This symmetry is broken by the formation of a chiral condensate. With maximum break- 
ing of chiral symmetry consistent with the Vafa-Witten theorem only an OSp(2\Nf + 2) 
subgroup remains unbroken. The Goldstone manifold is then given by the maximum Rie- 
mannian submanifold ]57]] of the coset space Gl(Nf + 2\2) /OSp(2\Nf + 2) with fermion- 



fermion integration domain given by U(Nf + 2)/0(Nf + 2) and boson-boson integration 
domain equal to U*{2, R) / Sp{2). 

An explicit representation of this manifold is given by 

U = V T G A V, (29) 

where 



G, 



' liVj+2 

1 j (30) 

V o 




and an explicit parameterization of V is given by 



V 



I UfF Oil Oi2 

fa a b |. (31) 

\ fa c d 



where Upp is an (Nf + 2) x (Nf + 2) unitary matrix, oli and fa are Grassmann variables 
and a, b, c and d are real numbers. 



4. The Scalar Susceptibility in ChPT 

In this section we compute the scalar susceptibility according to standard ChPT with 
Nf sea quarks and no valence quarks. The scalar susceptibility is defined by 

K ah = - f d'xd i yTTt a {x\-^—\y)t\y\-^—\x). (32) 
J ip + m Jp + m 

For equal quark masses, the Dirac operator is flavor diagonal. Replacing the spatial 
integrals by a sum over the Dirac spectrum, we obtain 



^=4?cd^- (33) 
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By introducing the isovector scalar source term s a qt a q the scalar susceptibility can 
be expressed in terms of a derivative of the QCD partition function (notice that the 
disconnected contribution vanishes in this case) 



K ab = d Sa d Sb hgZ. 



(34) 



Alternatively, the low-energy limit of the scalar susceptibility can be obtained from 
the effective chiral partition function. To lowest order in the momenta and the quark 
masses this partition function is determined uniquely by the pattern of chiral symmetry 



breaking |6T], p2| 



72 



(35) 



-S-lrid^Ud^U-^-^TTiMpiU+U- 1 )) 

dUe " L J 

lUeG/H 

The mass matrix also includes an isovector scalar source term. In all three cases, the 
matrix U is parametrized as 

U = e ^ Ta ' F ^ (36) 

where the sum is over the generators T a of the coset (notice that in this representation 
the mass matrix is diagonal for all three cases). The pion decay constant Fp is chosen 
such that we have the standard Gell-Mann-Oakes-Renner relation for all three cases with 
the convention TrT a T fe = 5 ab /2. For the three cases under consideration the appropriate 
parameters are given in Table 1. 



Table 1: The mass matrix and the pion decay constant for the three different patterns of chiral 
symmetry breaking. 



G/H 


M/3 






SU(N f ) 


ml Nf + s a t a 


F 


2 


SU(2N f )/Sp(2N f ) 


( (ml Nf + s a t a )/2 \ 
V (ml Nf + s a t aT )/2 ) 


F/V2 


1 


SU(N f )/0(N f ) 


(ml Nf + s a t a )/2 


F/V2 


4 



To second order in the 7r fields we find 



hu + U~ x ) = \-l S T, m T k T l . 
I bp 



(37) 



By differentiating the partition function with respect to the scalar source terms and 
calculating the loop integral we obtain 

4E 2 



4> r 
K ab = — Tr(t a T k T l )Tr(t b T m T n )(j rf 4 ^V fc (x)7r,(a:)7r m (y)7r n (y)) 1 _ loop 



F 4 
16tt 2 F 4 



Tr(t a {T\ T l })Ti(t b {T k } T 1 }) log(A/m). 



(3? 
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The only nontrivial task is to calculate the sum 

QW = ^Tr(t a {T k , T l })Tr(t b {T k , T 1 }) (39) 

Let us first consider the case (3 = 2. We already have seen that K ah is diagonal in isospin 
space. For the generators t a and t b we conveniently choose 

t a = t b = ^(5 tl 5 n -5 t2 5 j2 ). (40) 

The remaining diagonal generators are chosen as (consistent with the relation Tr t k t l = 

diag(^) = (V^l, -n,0 1 ^ 1: 0)/^2(n 2 + n), (41) 

n Nf—n—l 

whereas the off-diagonal generators are given by 

t k A = i(SikSji - S jk 5 u )/2, (42) 

and 

t k s = (SikSji + 8j k 5u)/2, (43) 

Notice that for (3 = 2 we do not distinguish between T a and t a . For the sum we thus 
obtain 

Q<» = I "f ' _i_ + ± {N , _ 2) = HOzl, (44) 

A^ 2 n(n + 1) 16 V ' ; 16 iV) ' V 1 

where the first term in this equation is due to the diagonal generators and the second 
term of this equation is due to the off-diagonal generators. 

As usual the case of (3 = 2 is the simplest case. The next simplest case is (3 = 4 with 
coset space given by SU(Nf)fO(Nf). In the sum ( f4"4"D the same diagonal generators enter 
as in the case of (3 = 2. However, only the symmetric off-diagonal generators contribute to 
this sum in this case, so that the contribution from the off-diagonal generators is exactly 
half of what it was for (3 = 2. This results in 



n=2 



n(n + 1) 



Finally we consider the case (3 = 1. In this case the generators of the coset have the 
structure |53J 

= (46) 
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where A is a Hermitian traceless matrix and B T = —B (both are Nf x Nf matrices). 
The sum over the generators is again calculated by distinguishing the diagonal generators 
from the off-diagonal generators. The diagonal generators are given by 

n = 4 ( ? !)> (47) 



V2 \ t 



D 



where the t a D are defined in (fil|). In order to satisfy the symmetry requirements, the 
off-diagonal generators have four nonzero entries in this case. Adding up all contributions 
we find 



where the first term is the contribution from the diagonal generators and the second term 
is the contribution from the off-diagonal generators. 
Our final result for K ab can be summarized as 

With the same assumption for the A and m dependence of the spectral density as in ||J] 
we find that the slope of the Dirac spectrum at the origin is given by 

p'(o) - s2 Wf-Wf + P) m 

MUj "l67r 2 F 4 (3N f • (5Uj 

The slope vanishes for Nf = 2 for all three values of (3. This is in agreement with instanton 
simulations [ETJ performed for the cases of (3 = 1 and (3 = 2. 



5. Partially Quenched Chiral Perturbation Theory 

The physics of QCD at low energies is completely determined by the spontaneous 
symmetry breaking pattern described above |)T], |63|. In general, the effective 



Lagrangian of partially quenched QCD is based on the spontaneous chiral symmetry 
breaking pattern of the graded Lie-groups G — > H. The Goldstone fields are given by 

U = exp(i2$/F /3 ), (51) 

where Fp is the pion decay constant and $ = <paT a with T a the generators of the maximum 
Riemannian sub manifold of G/H (notice that in this representation the mass matrix is 
diagonal for all three cases). We use the normalization StiT a T b = \5 ab . To lowest order 
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in the momenta and the quark masses, the Euclidean effective Lagrangian is obtained in 
the very same way as in Chiral Perturbation Theory pi , 32, 13, §]. This results in 



£ eff = 5str(^[/^C/- 1 ) - ^Str(^([/ + U- 1 )) + mg$g + ad,%d,%. (52) 

The last two terms in (|52|) represent the mass and kinetic terms of the super-?/ flavor- 
singlet field $o = Str($). They are introduced for the same symmetry reasons as in the 
three-color case |T3| , Below we are only interested in the case that m ~ ^qcd- in 
that case the term proportional to a is sub-leading, and it will be ignored from now on. 
The effective partition function is given by 

Z cS = [ dUe-I d4xC ^, (53) 

JU&G/H 

where the coset, the mass matrix and the pion decay constant for the three values of (3 
are given in the table below. 



Table 2: The mass matrix and the pion decay constant for the three different patterns of chiral 
symmetry breaking for the partially quenched effective partition function. 



G/H 



M, 



Gl(Nf + 1|1) 



Gl(2N f + 2|2) /OSp{2N f + 2|2) 



Gl(Nf + 2\2) /OSp(2\N f + 2) 




F/V2 



F/V2 



The mass matrix in the fermion-fermion block of the mass matrix Ain is given by 



M 2 


= diag 




• • • 


m s ,m v + J), 


Ma 


= diag 


(m s 


• • • 


m 8 , m v + J), 


M 4 


= diag 




Nf 

• • • 


m s , m v + J, m v + J) 








Nf 





(54) 
(55) 
(56) 



For the zero momentum component of the flavor singlet mass term, the same physical 
interpretation as in the three-color case holds ||: 



(57) 
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The factors in front of the first two terms in the effective Lagrangian ( p2|) have been 
chosen such that the Goldstone fields are normalized in the usual way. The masses of the 
Goldstone bosons are given by the usual Gell-Mann-Oakes-Renner relation 

M|. = {mi + rrij) S /F 2 . (58) 

with quark condensate in the chiral limit denoted by So- 

In the case of chiral perturbation theory for QCD it is possible to distinguish a domain 
where the partition function factorizes in a piece dominated by the zero momentum modes 
and a factor involving the nonzero momentum modes 0, |64|. The same is true for the 
partially quenched chiral Lagrangian. For valence quark masses below the Thouless energy 
given by §, || 1 



m c ~ — 2 , (59) 

the fluctuations of the zero modes in the the effective partition function dominate the 
fluctuations of the nonzero momentum modes. Then the integral over the zero momentum 
modes factorizes from the partition function and the calculation of H(m v ) is reduced to a 
group integral. This integral has been calculated explicitly for (3 = 2 with the result that 
Tj{m v ) is in complete agreement with chiral Random Matrix Theory. 



6. Calculation of T,(m v ) for (3 = 1 



From a computational point of view, two-color pqChPT is not very different from 
three-color pqChPT ||. To perform a calculation to one-loop order one has to construct 
the meson propagator. The effective Lagrangian discussed in previous section is based on 
the symmetric super-space Gl(2Nj + 2|2) /OSp(2Nf + 2|2). An explicit representation of 
$ in (51) in terms of a (2Nf + 4) x (2Nf + 4) matrix field is given by 



( 


i> 


Vv/2 


W 2 \ 






Xv/2 


X w /2 


Xv/2 


-Vv/2 


i(p v /V2 




V Xw/2 


-W2 


i(f> vw /2 


i(f> w /V2 J 



(60) 



where the (Nf+1) x (iVy + 1) ordinary matrices = $ and ip = —ip T contain the ordinary 
mesons made of quark and anti-quarks, the Nf + 1-vectors Xv, Xw, Vv an d f] w represent 
the fermionic mesons consisting of a ghost quark and an ordinary anti-quark, and finally 
the real fields (f> v , 4> w and 4> vw are the mesons with two ghost quarks. With the proper 
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normalization of the fields the Hermitian matrix field is given by 

/ 1)1 /2 



\w 2 ^ ••• 

and the anti-symmetric matrix ip is given by 

/ M2V2 

-l/>l, 2 /2y/2 



liJV/+1 /2 v / 2 \ 

<j>N f +l,N f +l/2 J 



(61) 



^ hNf+1 /2V2 



{ -4> hNf+1 /2V2 



lpN f ,N f +l/2V2 

-TpN f ,N,+i/2y/2 



(62) 



To second order in the pion fields there is no mixing between the off-diagonal Goldstone 
modes. However, because of the term Str 2 $, the situation is more complicated for the 
diagonal Goldstone modes. The propagator of these modes is obtained by inverting the 
quadratic form in the effective Lagrangian obtained by expanding the U fields up to second 
order. The mixing matrix of the diagonal mesons is given by 

/ 



V 



2 


2 


-iy/2 


-iy/2 


2 


2 


-iy/2 


-iy/2 


iy/2 ■■ 


• -iy/2 


-1 


-1 


■■ 


■ -iy/2 


-1 


-1 



(63) 



The propagator of the diagonal mesons is thus given by 

G^(p 2 ) = (T + m 2 M^)-\ 

where T is diagonal with 

p 2 + M 2 S for 1 < % < N f , 

p 2 + M 2 V for N f + 1 < % < N f + 3, 



r = 
1 n 



(64) 



(65) 



The matrix T + can be inverted by expansion in a power series in T 1 M^ and using 
the relation that 



(r -1 M (1) ) 2 = 2N f (p 2 + Ml)- 1 ^ 1 M (1) . 



By re-summing the geometric series we obtain 
G<V) ' ' 



_ m 2l M (l)l 



r ( p 2 + mi)- 1 + 2N f m 2 u r r 



(66) 



(67) 
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Below, we only need this propagator at the origin for 1 < i < Nf + 1 which is given by 

1 2m\ (P 2 + Ml) 



G 



(i) 



V ^ 

v 



(6* 



_ p 2 + Ml (p2 + M 2 )2 (p 2 + M 2 + 2jV/m 2) 

The quenched limit is obtained in the limit of infinite sea quark masses at fixed values 
of mo and m v . Standard ChPT power counting rules are recovered when the super- rj 
decouples, that is when the singlet mass mo is sent to infinity. It is easy to check that in 



this case the propagator (|68D in the sea quark sector is just the standard propagator in 
a quark basis with the singlet channel projected out. Finally, if both the valence and sea 
quark masses are set equal, we recover the standard result for two-color ChPT. 

The valence quark mass dependence of the chiral condensate is computed as in the 
three color case [jR], [!J. Since we are mainly interested in the non-analytic valence quark 
mass dependence of the chiral condensate, we will perform a one-loop computation, and 
we will ignore the analytic corrections coming from the 0(p 4 ) effective Lagrangian. 

The usual tadpole diagrams have to be evaluated. At finite volume, the pion propa- 
gator at the origin is 



(69) 



where the sum is over the momenta in the box. In the thermodynamic limit the sum can 
be replaced by an integral resulting in (in four Euclidean dimensions) 



A(M>) = 



(70) 



where A is the cutoff. 

It is now straightforward to calculate the valence quark mass dependence of the chiral 
condensate. Collecting the different contributions we find 



s (1) K) 



S 



— {2JV)A«) + G« - 2A(M 2 VV )} 



1 - 



N f F* 



2N 2 f A(M 2 vs ) - (N f + 1)A(M, 



2 > 
vv i 



+ {M 2 SS - Ml)d M ,A{Ml 



(71) 



The first term in the first equation is due to the off-diagonal mesons consisting out of a 
valence quark and a sea quark, the second term is due to the diagonal meson consisting 
out of two valence quarks, and the last term is due to mesons consisting out of a bosonic 
and a fermionic quarks. As usual there is an additional minus sign for fermionic meson 
loops. 
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7. Calculation of E(m v ) for (3 = 4 



For fermions in the adjoint representation the Goldstone manifold is determined by 
the pattern of chiral symmetry breaking according to 



Gl(N f + 2|2) -> 0^(21 JV/ + 2). 



(72) 



The Goldstone manifold is given by the maximum Riemannian submanifold of Gl(Nf + 
2|2)/0if>p(2|iV/+2). The effective Lagrangian has exactly the same structure as for (3 = 1, 
but in this case an explicit representation of the Goldstone modes is given by 



/ 



Aii/V2 



^1,^+2/2 
A/2 
V -ai/2 



A hNf+2 /2 a a /2 A/2 ^ 
^+2,^+2/^2 a Nf+2 /2 (3 Nf+2 /2 



^Nf+2/2 
-a Nf+2 /2 



#i/2 




(73) 



#1/2 y 



Here, A is a real symmetric matrix, 0i is a real variable and the and A are Grassmann 
variables. 

In order to do a one-loop calculation we need an explicit expression for the pion 
propagator. Because of the singlet mass term the inversion of the quadratic form involving 
the diagonal mesons is nontrivial. The propagator of the diagonal mesons is given by 



G^(p 2 ) = (T + m 2 M^)- 1 , 
where the mixing matrix is given by 



(74) 



and T is diagonal with 



P 2 + Ml 
P 2 + M 2 VV 



1 •■■ 1 -iy/2\ 

1 ••• 1 -iy/2 
-is/2 ■■■ -2 



for 1 < i < N f , 

for N f + 1 < i < N f + 3. 



(75) 



(76) 



The inverse in (f74|) can be calculated simply by expanding the expression as a geometric 
series in T^M^. If we use the relation 



N 



1 



V + Ml 



(77) 
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the geometric series can be resummed as 



r 1 + (p 2 + MD^Nfml °r 



(78) 



Below we only need the first Nf + 2 diagonal elements of the propagator matrix at the 
origin. They are given by 



(4) 



V ^ 

v v 



m 2 (p 2 + Ml 



p 2 + Ml {p 2 + M 2 )2(p2 + M 2 + Nfm T^ ■ 



(79) 



The calculation of the valence quark mass dependence of the chiral condensate proceeds 
in exactly the same way as for (3=1. As final answer we find 



s (4) K) = s c 



2F 2 



{iV / A(M s 2 J+2Glt ) -A(M„ 2 j} 



(80) 



8. Valence Quark Mass Dependence of the Chiral Condensate and the Dyson 
Index 



Before studying the valence quark mass dependence of the chiral condensate in five 
different limiting cases we write S(m B ) as a single formula valid for all three different 
values of (3. 

The result for the valence quark mass dependence of the chiral condensate for gauge 
fields in the fundamental representation and three or more colors was obtained previously 
T0| , H . In our present notation the result can be rewritten as 



m, 



1 - j- 2 {NfA(Ml) + Gg> - A(ilC)} 



(81) 



where G$ = G^j. With Dyson index (3 of the Dirac operator equal to /3 = 1 for QCD 
with two colors and fermions in the fundamental representation, (3 = 2 for three or more 
colors and fermions in the fundamental representation, and (3 = 4 for two or more colors 
in the adjoint representation, the result for the three different cases considered above can 
be written as 



So 



1 

F 2 



2Nf MMl) + G^-^A(Ml 



P 



(82) 



The diagonal propagator at the origin is given by 



1 

v 



E 



;i + cWK(p 2 + <) 



f + Ml {f + Ml){f + Ml + N f (l + 5^)ml) 



(83) 
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Notice that this propagator becomes independent of (3 in the limit mo — > oo when the 
usual ChPT counting rules prevail. Previously, a similar unified description was obtained 
for sum-rules for the inverse eigenvalues of the Dirac operator [[U|, [56|]. Below we will 
discuss five limiting cases for E(m„) obtained in (82). 

i) In the limit m v = m s <C uiq we should recover the mass dependence of the chiral 
condensate in standard chiral perturbation theory. In this case with M vv = M ss = M, 
the diagonal propagator reduces to 



N f 



N f 



-A(M 2 ), 



and the expression for E(m„) can be written as 



m, 



So 



1 - 



1 (N f -l)(l + 2N f /i3) 



F 2 



N f 



A(M 2 



(85) 



We observe that for Nf = 1 the result is independent of the Dyson index (3. Indeed, in 
this case the only axial symmetry is broken explicitly by the anomaly for all three values 
of (3. The factor (Nf — 1)(1 + 2Nf/(3) exactly counts the number of Goldstone bosons 
which is equal to the number of generators of SU(2Nf) / Sp(2Nf) for (3 = 1, the number 
of generators of SU(Nf) for (3 = 2, and the number of generators of SU(Nf)/0(Nf) for 
(3 = A. 

ii) Next we consider the case that both m v mo and m s mo- In this case the diagonal 
propagator is given by 



(86) 



This results in the valence quark mass dependence 



Ti(m r 



1 f 2Nf A / „ ,o \ ,N f 







In the thermodynamic limit, the singular part of £(M„) is then given by 



E(m, ; 



E f Nf 



%-k 2 N } F a \ (3 



{ o( m v + m s) lo f 



m v + m s 
2^ 



2 

+ (m s + (N f (l - -) -2)m v ) lo. 



m, 



where // = A 2 F 2 /2E . In the limit of m s — > this result simplifies to 



E(m v ) = E 



(89) 
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iii) In the range X m i n rn v <C 1/L 2 Aqcd <C m Q , and m s <C m t the propagator at the 
origin (|68| ) is dominated by the zero-momentum mode |)3j, [J5| |, [| resulting in 

1 



(90) 



M 2 V 

vv 

The valence quark mass dependence of the chiral condensate in the limit m s <C m v then 
reduces to 



,2N f 1 1 1 
■ (3 + 2 /Tr^VEo 



(91) 



iv) Next we consider the quenched limit which is obtained for m s ^> m v and m s ^> itlq. 
In this case the diagonal propagator at the origin is given by 

1 (1 + d/3,i)m 2 



G (p) = -y 



p 2 + M 2 v (p 2 + Ml 



(92) 



This results in the valence quark mass dependence 



S(m,; 



So 
So 



1 - — 2 {(1 - 2//3)A(M 2 J + (1 + 5 A1 )m^^A(M„ 2 j} 



1 



16tt 2 F 2 



(1 + 5 ft i)mo loj 



+ 0(771^ 



(93) 



where we have only quoted the infrared singular terms. We find that the chiral condensate 
diverges for m v — > 0. 

In the domain m v <C m s the expression (p8|) is identical to the quenched result (|93|) 
provided (1 + 5p^)mQ = M 2 s /Nf. This is in accordance with the Ward identity relating 
the sea quark mass and the topological susceptibility (u 2 ) = FjM 2 s V/2N f (l + 6p A ) §§, 
and the relation (|57|). 

v) In lattice QCD with staggered fermions the Lagrangian does not possess the axial 
£7(1) symmetry at finite lattice spacing (although it should be recovered in the continuum 
limit). A detailed discussion of quenched lattice simulations for the valence quark mass 
dependence of the chiral condensate in this case was given in [R7I] . The effective partition 
function is then given by ( |53"D with an effective Lagrangian without singlet terms. The 
diagonal propagator at the origin is thus given by 

1 



Gm ~V^p 2 + Ml 



(94) 



in the cases i), ii) and iv). In the quenched case, the term ~ m 2 , is absent and the 0(m v ) 
terms have to be taken into account. We then find 

:i-2//3)S . m v 
„ — m>v log — 



E(m„) 



1 - 



(95) 



8tt 2 F 4 

For (3 = 2 the coefficient of the chiral logarithm vanishes and higher order corrections 
become important |68fl . 
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9. Spectral Density in the Ergodic and Diffusive Domains 



In the ergodic regime, for valence quark masses m v <C F 2 /T, Q L 2 the microscopic 
spectral density is given by chRMT |7]], whereas for Aqcd 3> fUv \ m m ~ ^/L 4 T> the 
one- loop calculation of E(m v ) is valid. In this section we calculate the spectral densities 
in each of these domains. We first show that in the ergodic domain the one-loop result 
for E(m„) coincides with the chRMT result. 

The QCD partition function for two colors and fundamental fermions belongs to the 
chGOE (j3 = 1) universality class |TB[. In a sector of topological charge \v\ and in the 



chiral limit, the result reads [70 



Pt («) = 2[J2N f + \v\(u) - J2N f + \v\+l{u)J 2Nf+ \ v \- 1 {u)} 

1 f u 

+ 2 J 2Nf+W\( u )[ 1 ~J dtJ 2Nf +\ v \(t)]. (96) 

The QCD Dirac operator with adjoint fermions belongs to the chGSE ((3 = 4) univer- 
sality class. In this case the microscopic spectral density is given by J7I], |70 

pf GSE (u) = u[J 2 Nf+2H (2u) - J Nf+2H+1 (2u)J Nf+2Whl (2u)] 

1 f 2u 

(2u) / dtJ Nf+2H (t). (97) 

For completeness we also give the microscopic spectral density for the chGUE {(3 = 
2) JF2], ^1] which describes the Dirac operator of QCD with three or more colors and 
fundamental fermions, 

p? gu » = ^[4 /+ m(«)-^/+m+i(«)^ /+ h-i(«)]- 

(98) 

The valence quark mass dependence of the chiral condensate in terms of the micro- 
scopic variable x = m^VEo is obtained via the relation 

E(x) f°° , 2x , , \v\ , . 

where the last term is the contribution from the exactly zero eigenvalues of the Dirac 
operator [^] . In the case of the chGUE this integral can be perform analytically resulting 
in 

— — = x[I Nf+ \ v \(x)K Nf+ \ v \(x) + I Nf+ \ u \+i(x)K Nf+ \„\-i(x)]-\ 

Zjq X 

N f + \u\ \u\ N f 
= 1 J - Li + J_L + ... = i L + .... 100 

XX X 
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The asymptotic result for x ^> 1 for the chGOE and the chGSE can be easily derived using 
this result and the asymptotic result of the remaining terms. The leading 1/x correction 
can be obtained by expanding before integration. The only additional integral that is 
required is 



J^x) 



1. 



The asymptotic result in terms of the Dyson index is given by 



AN f + (3-2 
2px : 



(101) 



(102) 



which coincides with the result fl9T|) obtained from chiral perturbation theory. Notice, 
that the leading order asymptotic correction is independent of v. This is no longer true 
for higher order corrections and a summation over all v is required to compare with results 
obtained for fixed value of 9 (which is zero in our case) J73J . 

In the diffusive domain, for 1/L 2 Aqcd <C m v <C Aqcd the valence quark mass de- 
pendence of the chiral condensate is given by the one- loop result fl38|). The spectral 
density is given by the discontinuity of S(mt,) across the imaginary axis (p(X)/V = 



Disci 



£K)/2tt) 



As final result we then find for the spectral density of the 



Dirac operator with Nf flavors with quark masses m s <C A <C Aqcd 



pW 

V 



So 

7T 



1 + 



\Qtx 2 N 



F 4 j' 



(2N 2 f 



I A 







|A|Arctg^ - (N f (- - 1) + 2)tt|A| (103) 



N 2 



m 

A 2 + ml 



2m s log — 



The spectral density has a logarithmic infrared divergence when m s ^ 0. In the limit 
m s ^ A, the expression (|104j) reduces to: 



V 



7T 



1 + 



(N f -2)(N f + P)E 
16tt p N f F 4 



A I 



(104) 



The result for the slope for (3 = 2 was first obtained in JJJ. Remarkably, the slope vanishes 
for Nf = 2 for three values of (3. This is in agreement with instanton liquid simulations 
of the QCD Dirac spectrum [fil] . 

In the quenched case, the spectral density is found to be: 



P Q (X) 
V 



So 

TT 



(i + <WK , |A 



16tt 2 F 2 



log — 



(105) 



We find a logarithmic divergence of the spectral density for all three values of (3. This 
result is consistent with quenched instanton liquid and lattice QCD simulations ( 



74 



on liquid simulations and lattice QCD simulations |42 , 74| . For staggered fermions away 
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from the continuum limit this term has to be modified according to the discussion given 
in v) of the previous section. For the spectral density we then find 



10. Conclusions 

We have shown that the infrared limit of the QCD Dirac spectrum is completely 
determined by the pattern of chiral symmetry breaking. The spectral density follows 
from the low energy effective partition function of the partially quenched QCD partition 
function which in addition to the usual quarks contains valence quarks and their bosonic 
superpartners. In this paper we have analyzed the case of two colors with fundamental 
fermions and the case of fermions in the adjoint representation. In both cases the effective 
Lagrangian is based on a Riemannian super-manifold originating from chiral symmetry 
breaking to an ortho-symplectic graded Lie Group. As is the case for QCD with three 
or more colors with fundamental fermions the Goldstone manifold is characterized by a 
symbiosis of compact and noncompact degrees of freedom. 

Our one-loop calculation of the valence quark mass dependence of the chiral condensate 
for these two chiral symmetry breaking patterns completes an earlier calculation for QCD 
with three colors and fundamental fermions |J. Its discontinuity results in an explicit 
expression for the spectral density of the Dirac operator valid for A <C Aqcd- in all three 
cases we can distinguish two different domains in the Dirac spectrum which are separated 
by the Thouless energy. For valence quark masses below the Thouless energy the kinetic 
term decouples from the effective Lagrangian and the the one-loop result for the valence 
quark mass dependence of the chiral condensate is in complete agreement with chiral 
Random Matrix Theory. Agreement with chiral Random Matrix Theory to all orders has 
only been shown for (3 = 2 by the means of an exact calculation of the superintegrals. 
It would be interesting to perform this calculation for the mathematically much more 
complicated cases of (3 = 1 and (3 = 4 as well. For valence quark masses beyond the 
Thouless energy the kinetic term in the partially quenched effective partition function 
has to be taken into account. The results for the spectral density in this domain are 
consistent with the results for the scalar susceptibility obtained with standard ChPT. In 
particular we have obtained an analytical expression for the slope of the Dirac spectrum. 
Remarkably, the slope vanishes for two massless flavors in all three cases. 




(106) 



23 



Acknowledgements 

This work was partially supported by the US DOE grant DE-FG-88ER40388. One of 
us (D.T.) was supported in part by "Holderbank"-Stiftung. A. Altland, P. Damgaard, M. 
Gockeler, J. Osborn, P. Rakow, A. Smilga, M. Stephanov, T. Wettig and M. Zirnbauer 
are acknowledged for useful discussions. 



References 



[i] 

[2] 

[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 



T. Banks and A. Casher, Nucl. Phys. B169 (1980) 103. 

C. DeTar, Quark-gluon plasma in numerical simulations of QCD, in Quark gluon 
plasma 2, R. Hwa ed., World Scientific 1995. 

A.V. Smilga, Phys. Rep. 291, (1997) 1. 

A. Smilga and J. Stern, Phys. Lett. B318 (1993) 531. 

J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nuc. Phys. B540 (1999) 317. 



P.H. Damgaard, J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, |hep-th/9811212 
to be published in Nuc. Phys. B. 

J.J.M. Verbaarschot, Phys. Lett. B368 (1996) 137. 

A. Morel, J. Physique 48 (1987) 1111. 

C. Bernard and M. Golterman, Phys. Rev. D49 (1994) 486; C. Bernard and M. Goiter- 



man, 



hep-lat/931i07Tj 



[10] M.F.L. Golterman and K.C. Leung, |hep-lat/9711033| ; M.F.L. Golterman, Acta Phys. 
Polon. B25 (1994). 

[11] S. Chandrasekharan, Nucl. Phys. Proc. Suppl. 42 (1995) 475; S. Chandrasekharan 
and N. Christ, Nucl. Phys. Proc. Suppl. 42 (1996) 527; N. Christ, Nucl. Phys. B 
(Proc. Suppl.) 53 (1997) 253. 

[12] J.J.M. Verbaarschot, in Nonperturbative Approaches to Quantum Chromodynamics, 
D. Diakonov, ed., Gatchina 1995. 



[13] C. Bernard and M. Golterman, Phys. Rev. D46 (1992) 853; M. Golterman, |hep^ 



lat/9411005 , Chiral Perturbation Theory and the quenched approximation of QCD] S. 



Sharpe, Phys. Rev. D46 (1992) 3146. 



24 



[14 
[15 
[16 

[it; 

[18 

[19 
[20 
[21 
[22 

[23 
[24 
[25 
[26 

[27; 

[28 

[29 
[30 
[31 
[32 
[33 
[34 



S. Sharpe, Phys. Rev. D 56 (1997) 7052. 

E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. A560 (1993) 306. 

J. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531; Phys. Lett. B329 (1994) 351. 

E. Brezin, S. Hikami and A. Zee, Nucl. Phys. B464 (1996) 411. 

G. Akemann, P. Damgaard, U. Magnea and S. Nishigaki, Nucl. Phys. B 487[FS] 
(1997) 721. (1994) 487. 

A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Nucl. Phys. B479 (1996) 707. 
T. Guhr and T. Wettig, Nucl. Phys. B506 (1997) 589. 

A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Nucl. Phys. B506 (1997) 612. 

PH. Damgaard, Phys. Lett. B424 (1998) 322; B425 (1998) 151; G. Akemann and 
PH. Damgaard, Phys. Lett. B432 (1998) 390; Nuc. Phys. B528 (1998) 411. 

PH. Damgaard and S.M. Nishigaki, Phys. Rev. D57 (1998) 5299. 

J. Christiansen, |hep-th/9809194 . 



J.C. Osborn and J.J.M. Verbaarschot, Phys. Rev. Lett. 81 (1998) 268. 

M.A. Halasz and J.J.M. Verbaarschot, Phys. Rev. Lett. 74 (1995) 3920; M.A. Halasz, 
T. Kalkreuter and J.J.M. Verbaarschot, Nucl. Phys. Proc. Suppl. 53 (1997) 266. 

M.E. Berbenni-Bitsch, S. Meyer, A. Schafer, J.J.M. Verbaarschot and T. Wettig, 
Phys. Rev. Lett. 80 (1998) 1146. 

M.E. Berbenni-Bitsch, M. Gockeler, T. Guhr, A.D. Jackson, J.Z. Ma, S. Meyer, A. 
Schafer, H.A. Weidenmiiller, T. Wettig and T. Wilke, Phys. Lett. B438 (1998) 14. 

J.Z. Ma, T. Guhr and T. Wettig, Eur. Phys. J. A2 (1998) 87. 

T. Guhr, J.Z. Ma, S. Meyer and T. Wilke, Phys. Rev. D59 (1999) 054501. 

M.E. Berbenni-Bitsch, S. Meyer and T. Wettig, Phys. Rev. D58 (1998) 071502. 

K. Splittorff and A.D. Jackson, |hep-lat/9805UTg| ; K. Splittorf, |hep-th/98l0248 . 



R. Pullirsch, K. Rabitsch, T. Wettig and H. Markum, Phys. Lett. B427 (1998) 119. 

P.H. Damgaard, U.M. Heller, A. Krasnitz and T. Madsen, Phys. Lett. B440 (1998) 
129. 

25 



[35 
[36 

[37; 

[38 
[39 
[40 
[41 
[42 
[43 
[44; 
[45 

[46; 
[47; 

[48 
[49 

[50 
[51 



P.H. Damgaard, U.M. Heller and A. Krasnitz, hep-lat/981006. 

M. Gockeler, H. Hehl, P.E.L. Rakow, A. Schafer and T. Wettig, Phys. Rev. D59 
(1999) 094503. 



B.A. Berg, H. Markum, R. Pullirsch, |hep-lat/98120T0 



F. Farchioni, I. Hip, C.B. Lang and M. Wohlgenannt, |hep-lat/9812018 



R.G. Edwards, U.M. Heller and R. Narayanan, |hep-lat/9902TJ2l 



R.G. Edwards, U.M. Heller, J. Kiskis and R. Narayanan, |hep-th/9902TTT| . 

J.J.M. Verbaarschot, Nucl. Phys. B427 (1994) 534. 

J.C. Osborn and J.J.M. Verbaarschot, Nucl. Phys. B525 (1998) 738. 



A. Altland and B.D. Simons, |cond-mat/981113~4 



K. Takahashi and S. Iida, [hep-th/9903119 . 



B.L. Altshuler, I.Kh. Zharekeshev, S.A. Kotochigova and B.I. Shklovskii, Zh. Eksp. 
Teor. Fiz. 94 (1988) 343. 

H.A. Weidenmuller, Nucl. Phys. A 518 (1990) 1. 

K. Johnson, L. Lellouch and J. Polonyi, Nucl. Phys. B367 (1991) 675. 

R. Janik, G. Papp, M. Nowak and I.Zahed, Phys. Rev. Lett. 81 (1999) 264. 

K.B. Efetov, Adv. Phys. 32 (1983) 53; K.B. Efetov, Super symmetry in disorder and 
chaos, Cambridge University Press, (1997). 

T. Guhr, A. Muller-Groeling and H.A. Weidenmuller, Phys. Rept. 299 (1998) 189. 
G. Montambaux, in Quantum Fluctuations, Les Houches, Session LXIII, E. Gi- 



acobino, S. Reynaud and J. Zinn- Justin, eds., Elsevier Science, 1995, |cond- 
mat/960207j| . 



[52] F. Wegner, Z. Phys. B 49 (1983) 297. 

[53] J. Verbaarschot, H. Weidenmuller, and M. Zirnbauer, Phys. Rep. 129 (1985) 367. 
[54] M.E. Peskin, Nucl. Phys. B175 (1980) 197. 

[55] S. Dimopoulos, Nucl. Phys. B168 (1980) 69; M. Vysotskii, Y. Kogan and M. Shifman, 
Sov. J. Nucl. Phys. 42 (1985) 318; D.I. Diakonov and V.Yu. Petrov, Lecture notes in 
physics, 417, Springer 1993. 

26 



[56 

[5?; 

[58 

[59 
[60 
[61 

[62 
[63 

[64 
[65 
[66 
[67 

[68 
[69 
[70 
[71 
[72 
[73 
[74; 



A. Smilga and J.J.M. Verbaarschot, Phys. Rev. D51 (1995) 829. 
M. Zirnbauer, J. Math. Phys. 37 (1996) 4986. 

J.J.M. Verbaarschot, Lectures given at APCTP - RCNP Joint International School 
on Physics of Hadrons and QCD, Osaka, Japan, 1998 and the 1998 YITP Workshop 
on QCD and Hadron Physics, Kyoto, Japan, 1998, [hep-ph /9902394 . 



C. Vafa and E. Witten, Nucl. Phys. B234 (1984) 173. 

M.A. Halasz and J.J.M. Verbaarschot, Phys. Rev. D52 (1995) 2563. 

J. Gasser and H. Leutwyler, Ann. Phys. 158, 142 (1984); J. Gasser and H. Leutwyler, 
Nucl. Phys. B250, 465 (1985). 

H. Leutwyler, Ann. Phys. 235 (1994) 165. 

J. Gasser and H. Leutwyler, Phys. Lett. 188B(1987) 477; Nucl. Phys. B307 (1988) 
763. 

H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 5607. 

P. Hasenfratz and H. Leutwyler, Nucl. Phys. B343 (1990) 241. 

M. Shifman, A. Vainshtein and V. Zakharov, Nucl. Phys. B166 (1980) 493. 

M.E. Berbenni-Bitsch, M. Gockeler, H. Hehl, S. Meyer, P.E.L. Rakow, A. Schafer 
and T. Wettig, |hep-lat/99010T3 . 



P.E.L. Rakow, private communication. 

J.J.M. Verbaarschot, Nuc. Phys. B426 [FS] (1994) 559. 

T. Wettig, Habilitationsschrift Universitat Heidelberg (1998). 

T. Nagao and P.J. Forrester, Nucl. Phys. B435 (1995) 401. 

J. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70, 3852 (1993). 

PH. Damgaard, frep-th/ 9903096 . 



U. Sharan and M. Teper, [hep-lat/9808"0T7| ; frep-lat/ 98 12009 . 



27 



